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ABSTRACT 

O 

o 

(N ■ We study the probability that three GRBs (980425, 030329, 031203) are found 

within z = 0.17, given the luminosity functions consistent with the log — log 5* 
^yj • relationship for classical cosmological bursts {i.e., those observed by BATSE). 

Q> ' We show that, in order for the probability of these three low-z events to be 

^ ■ non-negligible (thus making it more likely that they belong to the same class 

of the classical cosmological bursts), the bursts' luminosity function must be a 
broken powerlaw. By reasoning in analogy with beamed AGNs, we show that 
^ observations are consistent with the expectations if GRB 980425 and GRB 031203 

Q\ • are indeed normal bursts seen sideways. Within this model, no bright burst 



> 

in 



o 

O ■ ~ 20 yr 



o 



within z = 0.17 should be observed by a HETE-like instrument within the next 



Subject headings: gamma rays: bursts 



1. Introduction 

After the discovery of the association of SN 2003dh with GRB 030329 (Stanek et ai, 
2003; Hjorth et al, 2003), and of SN 1998bw with 980425, it has been widely accepted that 
classical, cosmological GRBs arise from the simultaneous collapse of a SN to form a collapsar 
(Woosley 1994). Further evidence comes from the third nearby burst, GRB 031203, which 
also has been seen to be closely associated with a type Ic SN (Tagliaferri et al, 2004; Malesani 
et al, 2004). This model appears plausible for all long GRBs, even the distant ones, in the 
light of the similarity between these two objects and GRB 011121/SN2001ke (Bloom et al.. 
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2002; Garnavich et al. 2003), an object located at 2; ~ 0.36. Still, the fact that only some 
bursts' afterglows display the re-bumps now associated with the emergence of the optical 
contribution due to the underlying SN beckons the question of exactly which fraction of all 
GRBs (within the dctcctability range of SN rebumps, say z < 1) arc in fact associated with 
simultaneous SNc. In order to tackle this problem, we study whether the detection of GRB 
980425, 030329 and 031203, the only ones unequivocally associated with a Type Ic SN of 
extreme properties within a very small distance (z=0.17) from us, is a statistical anomaly, 
for the luminosity functions that are consistent with the GRB logN-logS distribution. 



2. Probabilities of detecting low- 2; events 

If the three low-rcdshift GRBs (980425, 030329 and 031203) arc really typical of the 
global GRB population, then their discovery within the current time and sky coverage, must 
be consistent with the local GRB explosion rate as deduced from the very large GRB samples 
made available by BATSE. A simple computation of the expected number of events within 
z — 0.17 (the redshift of GRB 030329) shows that this is not the case. The volume out to 
z = 0.17, in a = 0.7,0 = 1 cosmology with Ho = 71 km s"^ Mpc~^ is V ^ 1.4 Gpc^. 
Association of a GRB with a SN requires an accurate position, and thus only bursts revealed 
by the BeppoSAX WFC cameras, or HETE2 WXM telescopes, can be used. Inclusion of 
INTEGRAL would only make matters slightly worse. These two X-ray telescopes monitor 
Sb = 0.123 and Sh = 0.806 steradians respectively (Band, 2002). The total effective 
operation times for BeppoSAX and HETE2 are generously estimated as Tb = 4 yr and 
= 2 yr (L. Piro, private communication). If we call ho the observed local rate of GRBs, 
we find that the total expected number of GRBs located inside 2; = 0.17 is 

AT = noyM£±^ ^ 0.12 % (1) 

47r 0.5Gpc-^yr-i ^ ^ 

The probability of observing three bursts, assuming Poisson statistics, is then P3 = 2.7 x 

10-'Hno/(0.5 Gpc"^ JT-^)f. 

The major uncertainty in the above equation concerns ho, which takes on different 
values depending upon exactly which data property is fitted {V/Vmax, or the logA^ — logP 
relationship for BATSE data). For the local observed rate, Porciani and Madau (2001) 
consider three different star formation rates, and correspondingly find three values for ho, 
between 0.11 and 0.17 bursts Gpc~^ yr~^ (having converted their results to Ho = 71 km s~^ 
Mpc~^). Schmidt (2001) considers the same star formation rates, but fits different quantities, 
to obtain ho = 0.48, 0.51, 0.72 Gpc"^ yr-^ Perna, Frail and Sari (2003) find ho = 0.5 Gpc-=^ 
yr^^, while Guetta, Piran and Waxman (2003) find ho — 0.44 Gpc~^ yr~^. In the next 



-3- 



section, we will repeat the analysis of Guetta et aZ.(2003) with a more updated sample, 
to obtain rio = 1-1 Gpc~^ yr~^. With this value, the probability of a triple event would 
be P3 = 2.9 X 10~^(no/(l-l Gpc~^ yr~^))^. The hypothesis that the three nearby bursts, 
980425, 030329 and 031203, belong to the same group as the classical BATSE bursts, would 
thus be rejected with fairly high confidence. 

However, this is not correct. The local rate discussed above has been derived by all 
authors under the hypothesis that classical bursts exceed by far the luminosity of GRB 
980425 and they also exceed that of GRB 031203 (i.e. they considered a minimum luminosity 
for the LF L5 > 5 x 10^^ erg/sec, which from now on wc will assume to define the minimum 
luminosity for the classical population). Therefore, the derived probabilities are not self- 
consistent. An unified picture can only be obtained with a LF that includes all luminosities 
down to that of GRB 980425, and at the same time yields high probabihty of observing the 
three low- 2; events. This is what we will achieve in §4, after setting up the formalism and 
recomputing (for the more updated sample) the local rate for the classical population in §3. 



3. Luminosity function for classical GRBs 

We consider all 2204 GRBs of the GUSBAD catalog which contains all the long GRBs 
(Tgo > 2sec) (Kouvehotou et al., 1993), detected while the BATSE onboard trigger (Paciesas 
et al. 1999) was set for 5.5 a over background (resolution of 1024 ms) in at least two 
detectors, in the energy range 50-300 keV. We estimate Cmax/C'inin for each burst (where 
C'max is the count rate in the second brightest illuminated detector and Cmin is the minimum 
detectable rate) and find {V/V^^) = 0.335 ± 0.007. 

The luminosity function (LF) of the "classical" population of long-duration GRBs may 

be represented as a power law with lower and upper limits. Lb and L2, respectively. The 
local LF of GRB peak luminosities L, defined as the co-moving space density of GRBs in 
the interval log L to log L + d log L is: 

$„(L) = CoiL/L,y- , L^<L<L2 (2) 

where Co is a normalization constant. We stress that this LF is the "isotropic-equivalent" 
LF, it does not include the effects of beaming. 

We assume that GRBs trace the star-formation history, and adopt the recent SFR de- 
rived by Rowan- Robinson (1999; RR-SFR); this can be fitted with the expression Rgrb{z) — 
fio max(0.75, 0.75 z) 

The modeling procedure involves the derivation of the peak flux P(L,z) of a GRB of 
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peak luminosity L observed at redshift z: 

P(T .\ - ^ C{E,{l + z),E,{l + z)) 

where Dl{z) is the bolometric luminosity distance and C{Ei,E2) is the spectral energy 
distribution integrated between Ei = 50 keV and E2 = 300 keV. Schmidt (2001) finds that 
the median value of the spectral photon index in the 50-300keV band for long bursts GRBs 
is -1.6. We use this value for our analysis to include the k-correction. 

Objects with luminosity L observed by BATSE with a flux limit P^^ arc detectable to 
a maximum redshift Zmax(-^, -Piim) that can be derived from Eq. (3). We consider an average 
limiting flux = 0.25 ph cm^^s""*^ taken from the BATSE catalog. 



The number of bursts with a peak flux > P is given by: 

A^(> P) = J ^o{L)d\ogL 

^(^'^) Rgrb{z) dV{z) 



Jo 1 + z dz ^ ^ 

where the factor (1 + z)~^ accounts for cosmological time dilation and dV{z)/dz is the 
comoving volume element. 

To be consistent with previous calculations, we take the low luminosity cut-off Lj, ~ 
5 X 10^^ erg s~^, while the high luminosity, L2 ~ 5 x 10^^ erg s~^, is taken on the order of 
the maximum luminosity detected until the time of writing (Bloom et al. 2003). The slope 
a of the LP is constrained by fitting the model with the observed peak fiux distribution 
with a non-linear Levemberg-Marquart minimum method. To avoid problems with error 
correlation in a cumulative distribution like A^(> P) propagate in an unknown way, we use 
the differential distributions n(P) = dN/dP for this analysis. We find that the best fit is 
given for a ~ 0.72; the p- value of the fit is ~ 0.70 (see Fig.l upper panel). Both a and the 
normalization value are somewhat insensitive to the choice of below a value ~ 10^° erg 
s~^. This is mainly because GRBs with very low luminosity appear above the sensitivity 
limit of ~ 0.25 photons cm~^ s~^ in a very small volume around the observer. 

To obtain the observed local rate of GRBs per unit volume, no, we need to estimate 
the effective full-sky coverage of our GRB sample. The GUSBAD catalog represents 3.185 
years of BATSE full sky coverage implying a rate of 692 GRB per year. Using our LF we 
find ho ~ 1.1 Gpc~^yr~^. ^ 



^Note that for a different SFR like for example the SF2-SFR of Porciani and Madau (2001) we still find 
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4. Luminosity function for bursts seen from any direction and revised 

probabilities 



Wc arc thus left with the question of how to inchide GRB 980425 and GRB 031203 within 
our analysis. It is quite possible that these bursts belong to a different, local population of 
bursts, as occasionally claimed by some authors for GRB 980425 (Bloom et al. 1998). Still, 
there are at least two arguments suggesting that 980425 is just a normal burst seen sideways 
(Nakamura 1998; Eichler & Levinson 1999; Woosley, Eastman k Schmidt 1999). One is the 
exceptional similarity of the two underlying SNe in 980425 and 030329, within a class (that 
of extreme Type Ic SNe) known instead for its lack of common patterns. The other is that 
GRB 021211 (della Valle et al, 2003), and all bursts with the so-called optical rebumps are 
likely associated also to SN. 

For this reason, we consider a LF model which includes GRBs seen sideways. We follow 
the analogy with beamed AGNs, according to which (Urry and Shafer 1984) the total LF 
is simply a broken power-law, agreeing, at the bright end, with the LF of the objects seen 
face-on. We thus take: 



We impose here Li — -L980425, L2 — 5 x 10^^ erg s~^ and Lb to be the border value between 
GRBs seen sideways (weak bursts) and the ones seen face-on (classical bursts). We choose 
somewhat arbitrarily L;, ~ 5 x 10^^ erg s~^, and fit the same data as in the previous section 
with a and 7 as free parameters. Note that it has been proposed that X-ray flashes (XRFs) 
are simply GRBs viewed off-axis. If this interpretation is correct, then XRFs could be 
included in our analysis in the lower part of the LF broken power law (recalling that XRFs 
are generally underluminous compared to GRBs). The best-fit {a = 0.7, 7 = 0.1 and p- value 
~ 0.76) is shown in Fig. 1 (lower panel). ^ With this LF, the total local rate of events, down 
to the lowest luminosity (that of 980425) is now much higher, ho ~ 10 Gpc^'^yr^^, while 
the rate for L > Lf, remains unaffected. Using Eq(l), we find that the expected number of 
events within z — 0.17, and of luminosity even as low as GRB 980425, is Ni — 2.4, which is 
clearly highly compatible with the current observation of 3 events. 



a good fit and a local rate smaller by a factor ~ 2. We find a similar result also for the LF studied in the 
next session. 

^Note that the curve is almost identical to the one corresponding to the best fit LF obtained for the 
classical bursts shown in the upper panel. This is because to add a low luminosity (L < 5 x 10''^ erg/sec) 
tail to the LF does not aff'ect the logN-logS distribution (low luminosity GRBs appear above the BATSE 
sensitivity limit in a reduced volume as explained in section 4). 




Li < L < Lb 
Lb < L < L2 . 



(5) 
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A strong radio emission is expected from bursts seen sideways at ~ 1 yr delay (Waxman 
2004ab, Livio & Waxman 2000, see also Waxman 2004a,b and reference within) as the jet 
decelerates to sub-relativistic speed and its emission approaches isotropy. However no late- 
time rebrightening was detected for either GRB 980425 or GRB 031203 (Sodeberg, Frail 
& Wieringa 2004; Sodeberg et al. 2003, 2004). These results suggest that neither GRB 
980425 nor GRB 031203 were off-axis events, and instead were intrinsically sub-energetic 
GRBs. Therefore we also consider the possibility that both GRB 980425 and GRB 031203 
might have intrinsically low luminosity and all GRBs are seen on axis. The LF can still 
have an intrinsic break for some (unknown) reason, and therefore it can still be represented 
by eq.(5) (as in Schmidt 2001). In this case we can leave Lj, as a, free parameter in our fit 
and look for the best fit parameters a,'y,Lb. The results {a = 0.95, 7 = 0.4, Lb = 10^^ erg 
s~^ and p-value ~ 0.64) are shown in Fig. 1 lower panel. With this LF the total local rate 
of events is rio = 20 Gpc~^ yr"^. However, if all GRBs are seen on axis, then there is no 
natural explanation for the break in the LF. Indeed a single power law LF with Li = -L980425! 
L2 = 5 X 10^^ erg s^^ could be rapresentativc of these bursts. However in this case we get 
a very high local rate, rio ~ 100 Gpc~^ yr~^, which would predict about 24 events within 
the z — 0.17 volume (again using Eq.l). The probability of having only 3 occurrences is 
~ 10"^, which makes this scenario very unlikely. Therefore, we conclude that a broken 
powerlaw for the luminosity function is highly preferred to a single powerlaw, if the three 
low-^ events do belong to the same class as the classical ones. Also note that the rate for the 
single powerlaw would yield about 400 GEM (galactic events per million years), implying an 
implausibly large number of GRB remnants (Loeb & Perna 1998; Efremov et al. 1998) at 
any given time. 

Our model also allows us to compute the probability of finding one object (GRB 980425) 
as close as z = 0.008, which would be of course prohibitively small if we were to apply to this 
burst the statistics for classical GRBs discussed in §3. For our best-fit powerlaw model, we 
find that the probablity of bursts as luminous as GRB 980425 within the volume accessible 
to BeppoSAX and HETEII, is ^ 10"^, small but still not completely ruled out. For the 
detection threshold of Swift (Gehrels et al. 2004), this rate would be 0.05, while that of all 
bursts with L < L;, is ~ 2 for a 2-yr observation period. 
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2.8 I r 




Fig. 1. — Upper panel: The best-fit differential distribution of the classical GRBs vs. the 
observed one taken from the GUSBAD catalog. Lower panel: The best-fit differential dis- 
tribution of the GRBs in the case in which Lb marks the border between off-axis and on 
axis-GRBs (solid line) and in the case in which we leave Lf, as a free parameter (dashed line) 
vs. the observed one taken from the GUSBAD catalog. For both panels: The \ogio{CM/Cm) 
is divided into 14 bins, the first 11 of equal size and the remaining 3 with varying sizes 
chosen such that the number of bursts per bin in the observed sample is at least A^min = 40 
(in order to have reasonable statistics, so that the Poisson error will not be too large), and N 
is the observed {circles) or theoretical {lines; only the values at the center of each bin count) 
number of bursts in each bin. The edges of the bins are also plotted. 
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5. Conclusions 

We have shown in this Letter that the presence of three nearby GRBs does not pose a 
problem for the current view (according to which both distant and nearby bursts belong to 
the same class), if we modify the LF to include an extension to luminosities as low as that 
of GRB 980425. We did this by means of a minimum impact extension, i.e., by assuming, 
in complete analogy with beamed AGNs, that the LF of the bursts is a broken power-law, 
with the bright-end distribution equal to what we derived when we neglected the presence 
of faint, sideways bursts. As a by-product, we showed that this yields a non-negligible 
probability for the detection of GRB 980425-like event, which is also a new result. 

How could this picture change? Some authors have insisted that it is not possible to 
throw GRB 980425 in the same cage as classical GRBs, because of its unique properties in 
the radio band. Within our model, the local rate of nearby bright biirsts (L > 5 x 10"^^ erg 
s~^) observable by HETEII within z = 0.17 is UoV Sh/4:TT = 0.057 yr^^. This implies that 
we ought to observe the next such event ~ 20 yr from now. Should we see one significantly 
earlier than that then the argument suggesting a similarity between the distant, classical 
bursts and the nearby ones would have to be reassessed. 

We thank the anonymous referee for his/her very useful comments. We thank Da- 
vide Lazzati for his comments. DG acknowledges the NSF grant AST-0307502 for financial 
support. 
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